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A  cathode  half  cell  physics-based  model  for  a  St.  Jude  Medical  fabricated  silver  vanadium  oxide  (SVO) 
cathode  coin  cell  battery  was  constructed.  The  model  is  based  on  a  single  particle  Fields  second  law 
approach  with  the  open-circuit  potential  modeled  with  a  Redlich-Kister  equation.  By  assuming  that 
lithium  ions  intercalate  only  through  the  ends  of  the  cuboid  SVO  particles,  the  model  is  able  to  predict 
accurately  the  discharge  profile  of  experimental  cathode  half  cell  coin  cells. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Lithium  ion  silver  vanadium  oxide  (Li-SVO)  intercalation  elec¬ 
trode  batteries  are  one  of  the  most  commonly  used  cells  for 
modern  implantable  cardiac  defibrillator  devices.  The  SVO  reduc¬ 
tion  reaction  observed  during  discharge  has  been  characterized  by 
Gomadam  et  al.  [1  ]  as 

(*c  +yc  )Li+  +  (xG  +yc)e“  +Ag2+V45+On 

-  LiXG+yc+Ag2_Xc+VyG4+V4_yc5+011  +  XcAg° 

where  Xq  goes  from  zero  to  two  and  yc  goes  from  zero  to  four. 
Gomadam  et  al.  claim  this  reaction  can  be  treated  mathematically 
as  a  combination  of  two  parallel  reactions  [1  ]. 

xcLi+  +xGe-  +  Ag2+V45+On  -*  LiXc+Ag2_Xc+V45+0„  +xcAg° 

ycLi+  +yGe-  +  Ag2+V45+0„  LiyG+Ag2+VyG4+V4_yG5+0„ 

The  SVO  model  developed  by  Gomadam  et  al.  [1]  provided  an 
approach  to  modeling  of  Li-SVO  cathodes  by  using  Butler-Volmer 
equations  to  govern  electrode  kinetics  for  each  reaction  and  then 
linking  the  reaction  current  density  to  the  depth  of  discharge  (DoD) 
through  Faraday’s  law.  They  made  the  following  assumptions, 


1.  The  cathode  limits  battery  capacity,  such  that  there  is  excess  Li 
in  the  anode. 

2.  The  cathode  dominates  battery  resistance  and  contributions  to 
the  resistance  from  the  anode  and  separator  are  negligible. 

3.  The  cathode  is  kinetically  limited,  meaning  the  ohmic  and  mass 
transfer  resistances  are  negligible,  such  that  the  reaction  current 
is  uniform  throughout  the  porous  electrode. 

4.  The  effects  of  heat  generation,  degradation,  and  parasitic  reac¬ 
tions  are  negligible. 

5.  The  SVO  particles  are  cylindrical  in  shape. 

Their  model  was  limited  by  a  purely  empirical  fit  of  the 
open-circuit  potential  (OCP)  as  a  function  of  DoD  while  also  not 
addressing  potential  diffusion  limitations  from  mass  transfer  resis¬ 
tance.  Also,  they  indicated  that  in  their  hybrid  SVO-CF*  cathode 
system,  their  assumption  of  negligible  ohmic  resistance  was  valid 
up  until  a  particular  value  of  a  reaction  current  distribution  param¬ 
eter,  determined  by  a  ratio  of  the  effective  ohmic  resistance  to  the 
applied  current  density  and  based  on  the  cathode  design  param¬ 
eters.  A  later  paper  by  Gomadam  et  al.  [2]  then  expanded  their 
earlier  model  to  include  ohmic  resistance  effects.  The  paper  in  [2] 
does  not  mention  any  considerations  for  mass  transfer  limitations, 
however.  If  the  assumption  in  [1  ]  of  no  mass  transfer  resistances  is 
perpetuated  in  [2],  the  basis  for  their  limits  of  applicability  may  be 
incorrect. 

In  this  work,  a  proprietary  form  of  SVO  developed  by  St.  Jude 
Medical  (SJM)  with  a  one  reaction  system  is  used: 

7Li+  +  7e-  +  Ag2V45+0„  Li7+V24+V23+011  +  2Ag° 
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Nomenclature 

\ 

Redlich-Kister  coefficient  (J  mol-1 ) 

Ap 

particle  electrochemically  active  surface  area  (cm2) 

c 

dimensionless  lithium  ion  concentration 

cLi 

lithium  ion  concentration  (mol cm-3) 

cLi,0 

initial  lithium  ion  concentration  (mol cm-3) 

Ce 

electrolyte  phase  lithium  ion  concentration 
(mol  cm-3) 

^Lij.max 

maximum  intercalatable  lithium  ion  concentration 
(mol  cm-3) 

CLi  J,avg 

average  lithium  ion  concentration  (mol cm-3) 

^Li  J,surf 

surface  lithium  ion  concentration  (mol cm-3) 

Du 

concentration-dependent  diffusion  coefficient  of 
lithium  ions  into  SVO  (cm2  s-1 ) 

E 

electrode  potential  (V) 

F 

Faraday’s  constant  (96,485  C  mol-1  e-) 

Ge 

excess  Gibbs  free  energy  (J) 

H 

particle  half-height  (cm) 

Japp 

applied  current  (A) 

ij 

current  density  for  electrode  j  (A  cm-2) 

io 

exchange  current  density  (A cm-2) 

j 

positive  (p)  or  negative  (n)  electrode 

kj 

rate  constant  (cm2-5  mol-0-5  s-1 ) 

I< 

rate  terms  in  Eqs.  (23)-(25) 

L 

particle  length  from  center  to  one  end  (cm) 

M 

number  of  terms  in  Redlich-Kister  expansion 

ma 

active  mass  of  SVO  (g) 

MSv 

SVO  molecular  weight  (gmol-1 ) 

N 

number  of  nodes  in  method  of  lines 

Nli 

net  flux  density  of  lithium  ions  into  the  particle 
(molLi+  cm-2  s-1) 

ne 

number  of  electrons  transferred  during  total  reac¬ 
tion 

ideal  gas  constant  (8.314 J  mol-1  K-1 ) 

Ru 

rate  of  reaction  (mol  Li+  cm-3  s-1 ) 

Rohm 

ohmic  resistance  (Q) 

SU 

lithium  ion  stoichiometric  coefficient  for  total  reac¬ 
tion 

Sea 

electroactive  surface  area  (cm2) 

T 

temperature  (K) 

Tref 

reference  temperature  (K) 

t 

time  (s) 

t+ 

initial  time  (s) 

Un 

SVO  open-circuit  potential  (V) 

U° 

un 

standard  SVO  open-circuit  potential  (V) 

ULi 

mobility  of  Li+  (mol  Li+  cm2  J-1  s-1 ) 

Va 

volume  active  material  (cm3) 

Ve 

electrode  volume  (cm3) 

vP 

particle  electrochemically  active  volume  (cm3) 

V 

velocity  of  the  bulk  fluid  (cm  s-1 ) 

W 

particle  half-width  (cm) 

X-U,n,surf 

surface  depth  of  discharge  as  defined  by  Eq.  (19)  for 
the  negative  electrode 

*Li,  n,avg 

average  depth  of  discharge  as  defined  by  Eq.  (26)  for 
the  negative  electrode 

Z 

dimensionless  distance  from  center  of  particle 

Zu 

charge  number  of  the  reactive  species 

aaj 

anodic  transfer  coefficient  for  electrode  j 

acJ 

cathodic  transfer  coefficient  for  electrode; 

Pj 

symmetry  coefficient 

p 

bulk  SVO  density  (gem-3) 

T 

dimensionless  time 

Pj 

overpotential  of  electrode  j  (V) 

£ 

cathode  fraction  active  material 

Y Li,  n, surf 

activity  coefficient  of  intercalated  lithium 

<Pl 

solid-phase  potential  (V) 

cp2 

electrolyte-phase  potential  (V) 

$Li 

inflection  terms  in  Eqs.  (23)-(25) 

This  paper  extends  the  Gomadam  et  al.  model  by  using  the 
single  particle  framework  to  model  the  transport  within  the  parti¬ 
cle,  allowing  for  solid  phase  mass  transfer  limitations.  The  ohmic 
resistance  is  also  included  on  an  order  of  magnitude  basis  from 
experimental  SJM  data,  rather  than  assuming  a  negligible  ohmic 
resistance  within  a  range  of  applicability.  In  order  to  reduce  the 
dependence  on  empirical  fits,  the  model  presented  here  uses  the 
Redlich-Kister  (R-K)  equation  to  fit  the  open-circuit  profile  rather 
than  an  empirical  fit  as  performed  by  Gomadam  et  al. 

Atlung  et  al.  [3]  presented  a  single  particle  (SP)  model  for  an 
intercalation  electrode.  They  described  a  system  of  intercalation 
electrodes  using  an  electrolyte  salt  of  the  anode  material  and  an 
organic  solvent.  In  such  a  system,  cation  transport  was  affected  by 
diffusion  and  migration  effects,  where  the  diffusion  effects  domi¬ 
nated  in  the  solid  phase.  Cell  behavior  was  then  controlled  by  the 
electrolyte-electrode  interface  concentrations  of  the  ion.  With  this 
in  mind,  Atlung  et  al.  attempted  to  model  the  discharge  profile  of 
intercalation  electrode  batteries  based  on  the  cathode  particles  and 
their  geometry.  They  assumed  that  the  cation  and  electron  fluxes 
were  uniform  over  the  entire  particle  surface  and  that  all  cathode 
particles  were  in  contact  with  a  current  collector  matrix.  Under 
these  assumptions,  they  based  their  modeling  on  a  single  particle 
of  the  cathode  material  [3]. 

The  work  done  by  Atlung  et  al.  served  as  a  basis  for  the  SP  mod¬ 
eling  approach.  Haran  et  al.  [4]  expanded  on  the  work  by  Atlung 
et  al.  by  proposing  a  general  system  of  equations  for  solving  lithium 
intercalation  into  a  SP  intercalation  electrode.  The  work  done  by 
Haran  et  al.  included  work  by  Paxton  and  Newman  [5]  to  solve 
for  the  diffusion  coefficient  in  such  a  system,  but  such  work  is  not 
applied  in  this  paper  due  to  the  presence  of  plateaus  in  the  OCP  of  Li- 
SVO  that  Paxton  and  Newman  claimed  may  influence  their  method. 
This  paper  also  does  not  address  resistive  elements  expanded  on 
in  work  by  Zhang  et  al.  [6]  or  the  analytical  solution  for  diffusion 
applied  by  Guo  et  al.  [7]. 

Several  assumptions  are  made  here  when  using  the  SP  approach. 
As  presented  by  Santhanagopalan  et  al.  [8], 

1 .  The  concentration  of  lithium  ions  in  the  electrolyte  is  large  and 
independent  of  position  and  time. 

2.  The  ionic  conductivity  of  the  electrolyte  is  sufficient  enough  and 
the  flux  of  ions  small  enough  that  there  is  no  potential  drop  in 
the  solution  phase. 

3.  The  electronic  conductivity  of  the  porous  electrodes  is  large 
enough  and  the  current  small  enough  such  that  there  is  no  elec¬ 
tronic  potential  gradient  in  the  electrodes. 

4.  All  the  cathode  particles  are  uniformly  distributed  throughout 
the  cathode  and  are  of  the  same  geometry. 

Assumptions  1-4  guarantee  that  the  difference  in  the  potential 
of  the  electrolyte  and  active  material  is  independent  of  position. 
Thus,  the  flux  of  lithium  ions  into  each  particle  is  independent  of 
the  particle’s  position.  Therefore,  a  large  number  of  the  same  size 
particles  with  the  same  current  density  at  the  surface  of  the  particle 
can  be  used  to  represent  the  entire  cathode.  Also  assumed  for  this 
work: 
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Fig.  1.  SEM  image  of  pure  SVO  crystals  displaying  cuboid  geometry. 


5.  That  the  oxidation  of  lithium  metal  in  the  anode  occurs  close  to 

the  OCP  of  the  Li/Li+  reaction  and  that  there  is  a  large  excess  of 

lithium  metal  in  the  anode  compared  to  the  amount  of  lithium 

that  can  be  intercalated  into  the  cathode. 

The  main  advantage  of  applying  the  SP  method  over  other 
concentration-based  frameworks  is  the  SP  method  provides  a  sim¬ 
pler  system  of  equations  while  giving  accurate  results  for  the  C/100 
(253  p,A)  current  load  that  will  be  examined  in  this  work.  The  more 
significant  improvement  is  in  using  a  concentration-based  frame¬ 
work  instead  of  using  Faraday’s  law  to  link  the  reaction  current 
density  and  DoD  as  done  in  the  work  by  Gomadam  et  al.  Through 
the  concentration-based  framework,  one  can  include  mass  trans¬ 
fer  resistance  limitations  through  diffusion  into  the  particle,  not 
addressed  by  Gomadam  et  al.  This  offers  insight  into  how  the  con¬ 
centration  of  Li+  varies  inside  the  particle.  In  future  work,  this  can 
be  expanded  to  account  for  relaxation  from  pulse  effects  during 
pulse  discharge. 

2.  Modeling 

2.1.  Mode  of  intercalation 

A  SEM  image  of  SJM  SVO  particles,  Fig.  1,  indicates  that  the  par¬ 
ticles  are  cuboid  in  shape  instead  of  the  assumption  in  work  by 
Gomadam  et  al.  [1]  that  SVO  particles  were  cylindrical  in  geom¬ 
etry.  An  idealized  geometry  of  the  SJM  SVO  particle  is  shown  in 
Fig.  2  and  has  a  length  of  2 L,  a  height  of  2 H,  and  a  width  of  2  W.  SEM 
images  taken  by  Chen  et  al.  [9]  also  show  that  along  the  length  of 
these  particles  there  are  strata.  The  particles  themselves  have  been 
described  as  “tunnel-like”  [10].  The  work  of  Chen  et  al.,  combined 
with  observations  drawn  from  work  by  Takeuchi  et  al.  [11],  sup¬ 
ports  the  hypothesis  that  lithium  ions  intercalate  into  the  ends  of 
the  particle  along  the  length  of  these  strata.  However,  the  diffusion 
length  from  the  end  caps  to  the  center  of  the  particle  is  significantly 
larger  than  from  the  sides  of  the  particle  to  the  center,  which  may 
allow  lithium  ions  to  diffuse  through  the  sides  of  the  particle  as 
well. 

Thus,  two  limiting  cases  are  raised,  one  case  where  interca¬ 
lation  only  happens  from  the  ends  of  the  particle  and  one  case 
where  intercalation  only  happens  from  the  sides  of  the  particle. 
Mechanical  processes  used  by  SJM  to  grind  SVO  material  down  to 
smaller  particle  sizes  are  often  rated  by  their  ability  to  affect  par¬ 
ticle  length,  with  total  particle  length  varying  from  a  few  microns 
up  to  100  |jim.  Particle  width  and  height  also  generally  do  not  vary 


by  significant  margins,  typically  only  a  few  microns  in  size  and  are 
unaffected  by  the  grinding  process.  As  such,  in  the  interest  of  deter¬ 
mining  the  effect  of  mechanical  processing  on  discharge  behavior, 
this  work  will  concern  itself  with  the  limiting  case  where  interca¬ 
lation  only  happens  from  the  ends  of  the  particle.  Specifically,  the 
particle  will  be  modeled  such  that  Li+  intercalates  from  one  end  to 
the  center  of  the  particle.  The  particle  will  therefore  be  treated  as 
a  one-dimensional  slab  as  seen  in  the  simple  line  drawing  in  Fig.  3. 

2.2.  Fields  second  law 

The  material  balance  for  a  particle  is  described  by  Eqs.  (1)  and 

(2). 

y=-V'NLi+Ru  (1) 

Nu  =  -zuUuFcuV(fa  - DLi(cLi)VcLi  +  CLi V  (2) 

where  NLi  is  the  net  flux  density  of  lithium  ions  into  the  particle 
(mol  Li+  cm-2  s-1 ),  cL{  is  the  lithium  ion  concentration  in  the  particle 
(molLi+  cm-3),  RLi  is  the  rate  of  reaction  (molLi+cnrr3  s-1),  zLi  is 
the  charge  number  of  the  reactive  species,  uLi  is  the  mobility  of  Li+ 
(molLi+  cm2J-1  s-1),  tp\  is  the  solid-phase  potential  (J),  DLi  is  the 
diffusion  coefficient  of  lithium  ions  in  the  SVO  particle  (cm2  s-1 )  as 
a  function  of  lithium  ion  concentration,  and  v  is  the  velocity  of  the 
bulk  fluid  (cms-1).  Because  this  is  a  solid  phase,  v  =  0,  and  due  to 
assumption  3,  V</q  =0.  Also,  since  the  reaction  takes  place  at  the 
surface  of  the  particle,  ftLi  =  0.  This  reduces  the  material  balance  to 
Eq.  (3),  Fields  second  law. 

TF  -  -V-  <-Du(cu)Vcu)_DlJ(cu)  (^1  +  +  ^l) 

+  [V-DLi(cLi)]VcLi  (3) 


Eq.  (3)  can  be  simplified  by  assuming  that  diffusion  occurs  only 
along  the  direction  of  the  z-axis,  the  slab  approximation. 


9cu 

dt 


=  Du(Cu)  ■ 


92cLj 
9 z2 


9PLi(CLi) 

9cLi  ■  (9cLi/9z)2 


(4) 


At  some  initial  time,  t+,  the  concentration  of  lithium  ions  is  some 
initial  value,  cLi  0,  due  to  some  prior  discharge  to  test  for  possible 
cell  failure  before  beginning  normal  discharge  of  the  cell.  The  initial 
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Fig.  3.  Line  drawing  of  SVO  particle  studied.  Intercalation  from  one  end  to  center  only. 


concentration  is  left  as  an  adjustable  parameter  to  be  estimated 
using  nonlinear  least  squares. 


cLi(z,  t  =  4)  =  cLij0 


(5) 


Because  symmetry  was  assumed  there  is  no  flux  of  lithium  ions 
atz  =  0,  the  center  of  the  slab  [12]. 

%k0,t)  =  0  (6) 

At  the  surface  of  the  particle,  z  =  L,  the  flux  of  lithium  ions  into 
the  particle  is  proportional  to  the  current  density: 


dcLi  / ,  f  n  =  SU [app 
dz  ^  ;  neFDu(cLi)SEA 


(7) 


where  Iapp  is  the  applied  current  (A),  ne  is  the  number  of  electrons 
transferred  during  the  total  reaction,  sLi  is  the  stoichiometric  coef¬ 
ficient  of  Li+  in  the  overall  reaction  which  is  taken  to  be  seven  in 
this  paper  (sLi  =  7),  Fis  Faraday’s  constant  (96485  C  mol-1  e-),  and 
SEa  is  the  electroactive  surface  area  of  the  particle  (cm2). 

The  electroactive  surface  area,  SEA,  is  defined  as  the  electrochem- 
ically  active  surface  area  of  material  within  the  electrode,  or  the 
volume  of  active  material  within  the  electrode  (£eVe)  multiplied  by 
the  ratio  of  electrochemically  active  surface  area  of  the  particle  to 
particle  volume  (Ap/Vp)  [7]. 


Sea 


Ap£e  Ve 


(8) 


where  Ap  is  the  active  surface  area  of  the  particle,  Vp  is  the  volume  of 
the  particle,  Ve  is  the  volume  of  the  electrode,  and  ee  is  the  volume 
fraction  of  active  material.  The  volume  fraction  of  active  material 
is  defined  as  the  volume  of  active  material  within  the  electrode 
divided  by  the  total  volume  within  the  electrode.  If  the  mass  of 
active  material  and  the  bulk  density  of  the  material  are  known,  the 
volume  of  active  material,  Va,  can  be  determined  via  Eq.  (9). 


Eq.  (10)  defines  the  volume  fraction  of  active  material. 


o  _ 

£e“  Ve 


(10) 


Depending  on  particle  geometry,  the  surface  area  to  volume 
ratio  will  change.  For  slab  geometry,  where  the  particle  is  treated 
as  a  one-dimensional  entity  by  assuming  lithium  only  intercalates 
through  the  ends  of  the  particle,  and  examining  the  slab  from  its 
center  to  one  end  only,  Eq.  ( 1 1 )  applies. 


Ap  =  2(2  H  •  2  W)  VP=2L2W  ■  2  H 


Ap 

VP 


1 

L 


(11) 


Substituting  Eqs.  (9)-(  1 1 )  into  Eq.  (8)  expresses  the  electroactive 
surface  area  as  Eq.  (12). 


c  rna 

Sea=l7 


(12) 


Another  way  of  expressing  this  equation  is  that  the  electroac¬ 
tive  surface  area  is  the  mass  of  active  SVO  material  multiplied  by 
the  electroactive  specific  surface  area  of  the  section  of  the  particle 


under  examination.  If  one  were  to  assume  lithium  ions  intercalated 
into  the  SVO  particle  from  all  directions,  the  electroactive  surface 
area  to  volume  ratio  would  become  the  specific  surface  area  of  the 
entire  particle.  Since  the  endcap  intercalation  limiting  case  is  being 
investigated  in  this  paper,  only  the  specific  surface  area  of  the  ends 
of  the  particle  is  utilized  in  the  present  model. 

The  boundary  conditions  defined  by  Eqs.  (6)  and  (7)  with  the 
previously  defined  governing  equation,  Eq.  (4),  and  initial  con¬ 
dition,  Eq.  (5),  provide  the  basis  for  determining  lithium  ion 
concentration  within  the  particle  as  a  function  of  time  and  distance. 
In  order  to  increase  the  level  of  confidence  in  modeling  results,  the 
present  model  is  implemented  in  two  different  programs,  COMSOL 
Multiphysics  and  MATLAB.  COMSOL  allows  for  a  more  automated 
solution  of  systems  of  equations,  whereas  MATLAB  gives  greater 
freedom  in  programming  of  solution  steps.  In  order  to  allow  for 
changes  in  geometric  properties  in  COMSOL  without  having  to 
reconstruct  the  particle  geometry,  the  governing  equation,  initial 
condition,  and  boundary  conditions  can  be  expressed  in  dimen¬ 
sionless  terms  by  introducing  the  dimensionless  variables  in  Eq. 
(13)  [12]. 


C=^-l  z=v  r=Pu(cu)t 

cLi,0  L  L2 


(13) 


Applying  these  variables  to  the  governing  equation,  initial  con¬ 
dition,  and  boundary  conditions  gives  Eqs.  (14)-(17). 


C(Z,  r  =  r+)  =  0 

3C  _  32C  1  3DLi/3C\2 

dr  ~  dZ2  +DLi(cLi)  dC  \  dz) 

!<«•*)- 

r\ _ % IqppL 

dZ  ’  neFDu(cL[  )S£/\CLi0 


(14) 

(15) 

(16) 
(17) 


These  equations  can  then  be  placed  into  COMSOL  Multiphysics 
to  solve  for  the  dimensionless  concentration  as  a  function  of 
dimensionless  distance  and  dimensionless  time  after  specifying 
the  applied  current,  Iapp,  initial  concentration,  cLi0,  and  all  other 
parameters. 


2.3.  Reaction  kinetics 

The  generalized  electrochemical  reactions  for  intercalation  and 
deintercalation  at  a  solid/solution  phase  interface  as  described  by 
Guo  et  al.  [7]  are  written  as 

Li  —  6$  — >  Lr  +  e  +  0S 


Li  —  6>s  <—  Li+  +  e“  +  0S 

Intercalation 

where  Li -0S  is  an  intercalation  site  filled  by  a  lithium  ion  and  0S 
represents  an  unfilled  intercalation  site  [7].  The  rate  of  these  reac¬ 
tions  can  then  be  solved  with  Eq.  (18),  the  Butler-Volmer  equation 
(Appendix  A),  based  on  the  DoD  at  the  surface  of  the  particle,  xLij-surj, 
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defined  by  Eq.  (19)  and  the  electrode  overpotential,  rjj,  defined  by 
Eq.  (20)  and  neglecting  the  electrolyte  potential. 


h  ~  kj'^cLi,j,maxce'  xUj,Surf^  xli,j,surf) 


r 

"0.5  F 

0.5F 

x  <  exp 

wr,i 

-  exp 

W’1}_ 

} 

(18) 


v  _  ^Li ,j,surf 

xLi,j,surf  ~  ~  Oj,n,  max 

^Lij^max 


s_uP_ 

MSv 


Vj  ~  Ej  Uj 


(19) 

(20) 


where  subscript  j  refers  to  the  positive  (p)  or  negative  (n)  elec¬ 
trode,  ij  is  the  current  density  (A cm-2),  kj  is  the  rate  constant 
(cm2-5  mol-1  s-1 ),  Cnjjnax  is  the  maximum  concentration  of  lithium 
ions  that  can  be  intercalated  into  the  particle  (mol cm-3),  ce 
is  the  concentration  of  lithium  ions  in  the  electrolyte  phase 
( 1 E-03  mol  cm-3 ),  p  is  the  bulk  density  of  SVO  (4.789  g  cm-3  [  1  ] ), 
MSv  is  the  molecular  weight  of  SVO  (595.5  g  mol-1 ),  Ej  is  the  elec¬ 
trode  potential  (V),  and  Uj  is  the  OCP  (V)  referenced  to  a  Li/Li+ 
reference  electrode.  Because  a  Li/Li+  reference  electrode  is  used  for 
Uj,  the  anode  is  lithium  metal,  and  the  current  used  is  so  small, 
the  anode  potential  is  assumed  to  be  negligible.  Thus,  the  cathode 
potential  plus  the  ohmic  resistance  is  the  cell  potential  and  anodic 
terms  are  neglected.  Including  the  dimensionless  variables  from  Eq. 
(13)  and  the  ohmic  resistance,  Eq.  (18)  can  be  rewritten  as  Eq.  (21 ). 


En  —  lappRohm  +  EJn  +  q  ,-p  U  sinh 


2Iat 


S£ACe0-5[CLi,o(C(l,r)  +  l)]  •  [Cu.„fmax-Cu,o(C(l,T)  +  l)] 


(21) 


where  C(l,r)  is  the  dimensionless  concentration  at  the  surface  of 
the  particle  and  xu<n<smf=  cLi>0(C(l,  r)  + 1  )/cLi 

,n, max- 


2.4.  Open-circuit  potential  modeling 


The  OCP  of  Li-SVO  changes  over  the  course  of  DoD  such  that 
without  an  accurate  model  of  the  OCP,  the  discharge  profile  can¬ 
not  be  properly  modeled.  The  way  in  which  the  OCP  changes  with 
DoD  at  the  surface  of  the  particle,  however,  is  a  property  of  how 
the  material  was  fabricated,  requiring  fitting  of  an  OCP  model  to 
experimental  data.  In  literature,  a  purely  empirical  fit  to  the  data 
is  sometimes  used,  such  as  in  work  by  Gomadam  et  al.  [1].  The 
fit  in  [1]  is  peculiar  in  that  it  splits  the  OCP  data  curve  into  two 
separate  curves,  each  modeled  with  different  equations  that  are 
designed  to  be  run  simultaneously.  The  resulting  OCP  curve  from 
the  two  equations  working  together,  however,  does  not  match  up 
identically  with  the  source  data. 

In  order  to  minimize  the  level  of  empirical  analysis  and  match 
more  closely  with  available  OCP  data,  another  way  of  fitting  the  OCP 
is  used  in  this  work.  Karthikeyan  et  al.  [  1 3  ]  showed  that  the  OCP  of  a 
lithium  intercalation  electrode  can  be  accurately  modeled  using  the 
R-K  equation  to  estimate  the  excess  Gibbs  free  energy  (Appendix 
B).  The  resulting  OCP  equation  when  used  in  Eq.  (21)  is  given  by 
Eq.  (22),  where  DoD  is  defined  by  Eq.  (19)  as  the  Butler-Volmer 
equation  is  only  valid  at  the  surface  of  the  particle  [13]. 


FU„  =FU°  +RgTref  In 


1  Xn,n,surf 
xLi,n,surf 


M 


i  \k+l  ^xLi, n, surf  xLi,n,surf) 

iZxLi,n,surf  ~  1 ) 


i^xLi,n,surf  1) 


1  -k 


]} 


(22) 


The  advantage  of  using  the  R-K  fit  over  a  purely  empirical  fit  is 
the  R-K  equation  provides  physical  insight  for  why  the  OCP  curve 


Depth  of  discharge 

Fig.  4.  Comparison  of  the  open-circuit  potential  modeled  where  non-idealities  from 
activity  are  ignored  (Nernst  equation)  and  accounted  for  (R-I<  equation). 


changes  over  the  DoD,  relating  the  OCP  behavior  to  electrochemi¬ 
cal  phenomena.  Fig.  4  gives  a  comparison  of  the  OCP  of  the  SJM  cell 
under  study  where  the  non-idealities  due  to  activity  are  ignored, 
reducing  the  equation  to  the  classical  Nernst  equation,  and  where 
the  non-idealities  are  accounted  for  in  the  full  R-K  equation.  Work 
done  by  Colclasure  and  Kee  [14]  demonstrate  how  the  R-K  method 
could  be  expanded  to  estimate  the  activity  coefficient  and  then 
applied  that  to  determine  the  exchange  current  density. 

2.5.  Parameter  estimation 

The  R-K  coefficients  Ak  and  U°  are  considered  unknown  param¬ 
eters  in  the  model.  Since  SVO  is  sensitive  to  the  manner  in  which  it 
is  produced,  the  open-circuit  profile  is  distinct  for  each  method  of 
preparation.  Fig.  5  is  provided  as  an  example  of  how  the  OCP  profile 
for  SVO  can  differ  widely  based  on  the  manner  in  which  the  material 
was  produced  [1,15].  In  order  to  allow  the  R-K  equation  to  model 
the  OCP  based  on  input  data  from  a  desired  formulation  of  SVO,  the 
R-K  coefficients  are  considered  adjustable  parameters  and  deter¬ 
mined  by  performing  a  nonlinear  least  squares  fit  to  experimental 
OCP  data  using  MATLAB  and  the  MATLAB  routine  LSQNONLIN. 

Focusing  on  fitting  the  first  70%  DoD,  one  can  observe  from  Fig.  6 
that  there  is  a  difference  between  OCP  and  experimental  discharge 
data  that  is  larger  in  the  first  40%  of  discharge  in  comparison  to 
subsequent  regions.  This  has  been  reported  previously  within  lit¬ 
erature  and  attributed  by  Crespi  and  West  to  temporary  insertion 
of  lithium  ions  into  an  alternative  set  of  sites  due  to  the  slow  rate  of 
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Table  1 

Modeling  constants. 


Parameter 

Units 

Value 

Cu,n,max 

mol  Li+  cm-3 

5.629E-02 

Ce 

mol  Li+  cm-3 

IE-03 

lapp 

A 

-253E-06 

L 

cm 

25E-04 

ma 

g 

0.0803 

M 

- 

15 

Msv 

gmol-1 

595.5 

N 

- 

50 

ne 

mole-  mol-1  SVO 

7 

Sea 

cm2 

6.707 

T 

K 

310.15 

Tref 

K 

310.15 

ft 

- 

0.5 

P 

gem-3 

4.789 

5  Li 

mol  Li+  mol-1  SVO 

7 

Fig.  6.  Raw  SJM  OCP  and  C/100  (253  (jlA)  SJM  coin  cell  discharge  data. 


the  thermodynamically  preferred  reaction  [16].  In  this  paper,  the 
greater  difference  between  OCP  and  experimental  cell  potential  in 
the  first  40%  DoD  compared  to  the  rest  of  discharge  is  used  as  an 
aid  in  fitting  the  diffusion  coefficient,  DLi,  and  rate  constant,  kn.  The 
rate  constant  is  allowed  to  vary  over  the  extent  of  discharge.  The 
justification  behind  this  is  the  reduction  of  silver  to  silver  metal 
appears  to  have  a  significant  effect  on  cell  performance,  also  noted 
in  work  by  Gomadam  et  al.  [1  ].  Eq.  (23)  is  used  to  vary  the  rate  con¬ 
stant  for  this  work  and  was  derived  based  on  best  apparent  fit  to 
the  discharge  data. 


fen(Cu)  =  /<n,l  + 


1<n,2  ~  1<n,  1 


1  T  exp[  K]  ((Cy , n, surf  /  Cli, n, max)  ^Li,l)] 

kn,3  —  2 

1  +  CXp[— /^((Cy  n  slirf/CLi  n  max)  —  ^Li, 2 ) ] 


(23) 


where  the  I<  terms  define  how  quickly  the  equation  changes 
between  values  and  the  0Li  terms  refer  to  the  inflection  points 
in  the  resulting  curve.  Similarly,  the  diffusion  coefficient  DLi  has 
been  shown  in  literature  to  vary  over  the  course  of  discharge.  Data 
for  fitting  the  diffusion  coefficient,  however,  are  insufficient  within 
the  literature.  Results  from  Ramasamy  et  al.  [17]  show  a  diffusion 
coefficient  that  monotonically  decreases  with  DoD.  Results  from 
Takeuchi  et  al.  [18],  however,  report  a  diffusion  coefficient  that 
experiences  irregular  spikes  and  dips  as  discharge  proceeds.  The 
irregular  behavior  of  the  diffusion  coefficient  during  discharge  is 
also  supported  by  results  from  Lee  and  Popov  [19].  One  should 
note  that  the  results  obtained  by  the  Ramasamy,  Takeuchi,  and  Lee 
groups  do  not  agree.  This  would  imply  a  strong  dependence  on  the 
specific  mix  of  SVO  used.  Since  literature  sources  are  inconsistent 
and  there  are  insufficient  experimental  diffusion  coefficient  data 
available  to  develop  a  detailed  equation  for  the  SVO  mix  of  interest, 
an  empirical  equation  is  used.  Eq.  (24)  is  a  function  of  Li+  concen¬ 
tration  in  the  particle  and  was  developed  along  with  Eq.  (23)  from 
best  apparent  fit  to  discharge  data. 


^Li(cLi)  =  DU,\  +  Y 


DU,2  -  ^Li,l 


+  exp[-/C3((cLi(z)/cLi?n 

,n 

DLi,3  -  DU,2 


)  -  ^Li,3  )1 


1  +  exp[-/C4((cLi(z)/cLi!n?max)  -  0LiA)] 


(24) 


MATLAB  and  LSQNONLIN  can  be  used  to  fit  the  parameters  to  the 
data  with  a  nonlinear  least  squares  fit.  The  rate  constant,  diffusion 
coefficient,  and  initial  concentration  were  fitted  simultaneously. 
Eqs.  (4)-(7)  were  solved  using  the  method  of  lines  with  50  node 


points  in  MATLAB  (Appendix  C)  and  the  finite  element  technique 
in  COMSOL. 

The  assumption  is  made  in  work  by  Gomadam  et  al.  that  the 
effect  of  ohmic  resistance  is  negligible  at  low  current  loads.  In  this 
work,  the  ohmic  resistance  is  included  and  allowed  to  vary  over 
depth  of  discharge  through  Eq.  (25).  Since  resistance  data  were 
unavailable  for  the  fitted  coin  cell,  using  Eq.  (25)  with  the  parame¬ 
ters  defined  in  Table  3  allows  for  an  ohmic  resistance  on  the  order 
of  magnitude  seen  in  similar  SJM  SVO  coin  cells.  Additionally,  all 
values  of  I<  and  0Li  used  in  Eqs.  (23)-(25)  were  set  arbitrarily  to 
roughly  represent  observed  reaction  steps  in  work  by  Leising  et  al. 
[10]. 


Rohm(C  Li)  —  Rohm,  1  + 


K 


ohm, 2  Rohm,  1 


1  +  exp[  /^((Cy  ,n , surf /^Li  ,n  ,maX)  ^Li ,  5  )  ] 


+ 


Rohm, 3  Rohm, 2 

1  +  CXp[— JCg((Cy  n  sur//CLi?n,maX)  —  ^Li, 6 ) ] 


(25) 


3.  Results 

The  model  developed  in  this  work  was  applied  to  a  coin  cell 
described  in  Table  1  with  experimental  discharge  data  shown  in 
Fig.  6.  The  OCP  of  the  cell  is  independent  of  any  changes  to  SVO 
mass  or  particle  length  and  is  modeled  using  MATLAB  to  fit  the  R-K 
expansion  out  to  seventeen  terms,  including  U°  and  A0.  The  result¬ 
ing  fit  is  plotted  in  Fig.  7  and  the  R-K  coefficient  values  tabulated  in 
Table  2.  The  fit  extrapolates  to  the  final  reported  DoD  the  coin  cell 
was  modeled  to.  Part  of  the  attractiveness  in  using  the  R-K  equation 


Fig.  7.  Modeled  OCP  compared  to  experimental  SJM  OCP  data. 
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Table  2 

Redlich-Kister  coefficients  for  SJM  SVO  determined  through  parameter  estimation. 
U°  is  in  V  and  the  R-K  coefficient  terms  are  in  units  J  mol-1 . 


Parameter 

Value 

u„° 

2.608 

Ao 

— 4.576E+04 

Ai 

6.621E+03 

A2 

-4.578E+04 

A3 

-4.700E+04 

a4 

3.298E+05 

A5 

-1.872E+05 

A6 

-2.034E+06 

Ay 

5.935E+05 

As 

6.575E+06 

Ag 

5.883E+05 

A10 

— 1.097E+07 

An 

-3.529E+06 

A12 

8.993E+06 

A13 

4.067E+06 

Al4 

— 2.877E+06 

A15 

-1.521E+06 

0  0.2  0.4  0.6  0.8  1 

Depth  of  discharge  x  u,n,avg 


Fig.  8.  Activity  coefficient  determined  using  estimated  Redlich-Kister  parameters. 

is  that  it  can  be  related  back  to  the  activity  coefficient  (Appendix 
B)  [13,14].  The  activity  coefficient  is  plotted  versus  average  particle 
DoD,  defined  by  Eq.  (26),  in  Fig.  8. 


Lipquvg  «  v  7 

LLi,n,max 

In  work  done  by  Leising  et  al.  [10],  the  claim  is  made  that  elec¬ 
trochemical  activity  of  a  Li-SVO  cell  is  dominated  by  reduction  of 
silver  to  silver  metal,  which  lasts  until  around  34%  DoD.  This  is 
represented  in  Fig.  8  as  a  very  low  activity  coefficient,  meaning 
there  is  a  very  strong  tendency  for  Li+  to  intercalate  into  SVO.  In 
the  region  from  34%  to  54%  DoD,  Leising  et  al.  claims  the  electro¬ 
chemical  activity  switches  to  a  purely  V5+  to  V4+  reduction,  reflected 
as  a  higher  activity  coefficient  in  Fig.  8.  After  54%  DoD,  Leising  et  al. 
states  there  is  a  competition  between  reducing  V5+  to  V4+  and  V4+  to 
V3+.  The  OCP  drops  gradually  in  this  region  up  until  around  74%  DoD, 
when  a  more  dramatic  drop  in  potential  begins  to  occur,  leading  to 
the  conclusion  that  the  V4+  to  V3+  reduction  becomes  dominant  at 
this  point.  In  Fig.  8,  the  activity  coefficient  also  rapidly  begins  to 
approach  unity  after  around  74%  DoD.  Thus,  the  activity  coefficient 
is  demonstrating  the  sudden  drop  in  useful  capacity  beyond  74% 
DoD  as  lithium  ions  become  more  difficult  to  intercalate  into  SVO. 
These  results  provide  a  basis  for  using  the  74%  DoD  mark  as  an  end 
of  life  indicator  for  the  cell,  which  is  also  used  by  SJM. 

Both  MATLAB  and  COMSOL  can  now  be  used  to  show  the  concen¬ 
tration  and  discharge  profiles  with  both  time  and  position  within 
the  particle.  First  consider  the  previously  described  coin  cell  under 
constant  discharge.  MATLAB  is  used  to  perform  parameter  estima- 


Fig.  9.  Fitted  MATLAB  model  compared  to  SJM  OCP  data  and  C/100  (253  pA)  SJM 
coin  cell  data. 


Table  3 

Estimated  discharge  fitting  parameters  for  SJM  discharge  data. 


Parameter 

Units 

Value 

Qj.o 

mol  Li+  cm-3 

7.506E-05 

kn,  1 

cm2-5  mol-0-5  s-1 

3.781E-06 

kn,  2 

cm2-5  mol-0  5  s-1 

2.662E-07 

kn,  3 

cm2-5  mol-0-5  s-1 

1.873E— 06 

f\i,l 

cm2  s_1 

3.897E-10 

Du,2 

cm2  s_1 

7.969E-11 

Du,3 

cm2  s-1 

2.570E-08 

Rohm.l 

Q 

13.07 

Rohm,  2 

Q 

10.07 

Rohm, 3 

Q 

7.063 

K 1 

- 

50 

I<2 

- 

10 

I<3 

- 

50 

IQ 

- 

100 

K 5 

- 

100 

IQ 

- 

30 

0Li,l 

- 

0.100 

0\a,2 

- 

0.380 

#Li,3 

- 

0.090 

#Li,4 

- 

0.340 

#Li,5 

- 

0.150 

#Li,6 

- 

0.300 

tion  and  fit  the  model  to  the  experimental  data.  The  resulting  fit 
is  shown  in  Fig.  9,  where  DoD  is  based  on  Eq.  (26),  with  param¬ 
eter  values  in  Table  3.  The  modeled  fit  very  closely  matches  the 
experimental  data  except  for  the  region  from  5%  to  10%  DoD.  The 
discrepancy  between  the  model  predictions  and  the  experimental 
data  is  attributed  to  inaccuracies  in  how  the  rate  constant,  plot¬ 
ted  in  Fig.  10,  and  diffusion  coefficient,  plotted  in  Fig.  11  on  a  log 
scale,  are  determined.  Fig.  12  plots  the  diffusion  coefficient  calcu¬ 
lated  from  Eq.  (24)  on  a  log  scale  versus  cell  potential  and  compares 
those  results  to  results  reported  by  Lee  and  Popov  [19].  One  should 
note  that  the  SVO  used  in  the  work  by  Lee  and  Popov  may  be  con¬ 
siderably  different  from  the  SJM  SVO  studied  in  this  paper,  and  also 
approaches  the  diffusivity  from  an  electrode  basis  rather  than  a 
single  particle  that  has  been  limited  to  intercalation  from  one  end 
only.  From  the  comparison,  a  conclusion  could  be  drawn  that  the 
diffusion  coefficient  as  calculated  by  Eq.  (24)  should  start  at  a  value 
several  orders  of  magnitude  higher  and  then  decrease  as  the  DoD 
approaches  10%.  Increasing  the  complexity  of  the  empirical  equa¬ 
tions  utilized  is  expected  to  resolve  the  discrepancy.  Fig.  13  shows 
how  the  concentration  within  the  particle  varies  across  the  length 
of  the  particle  at  a  single  point  in  time.  This  leads  to  the  conclu¬ 
sion  that  on  the  chosen  current  load,  the  concentration  across  the 
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Fig.  10.  Rate  constant  varying  with  discharge  based  on  Eq.  (23)  and  estimated 
parameters. 


Fig.  11.  Diffusion  coefficient  varying  with  discharge  based  on  Eq.  (24)  and  estimated 
parameters. 


Fig.  12.  Diffusion  coefficient  as  calculated  from  Eq.  (24)  compared  to  results  from 
Lee  and  Popov  [19]. 


Dimensionless  distance  from  center  Z 

Fig.  13.  Concentration  profile  in  a  C/100  (253  |xA)  coin  cell  particle  with  changing 
distance  from  center. 


Fig.  14.  Comparison  between  Gomadam  et  al.  and  MATLAB  fits  to  Gomadam  et  al. 
1000  |xA  discharge  data  [1]. 

length  of  the  particle  could  be  considered  constant.  Even  at  the 
low  253  |jiA  current  load  chosen,  however,  the  effect  of  the  con¬ 
centration  gradient  across  the  particle  on  the  discharge  profile  and 
estimated  parameter  values  was  determined  to  be  significant. 

In  order  to  be  confident  in  the  validity  of  the  numerical  results, 
the  MATLAB  model  results  were  also  compared  against  results 
generated  by  COMSOL  Multiphysics,  which  numerically  solves  for 
Fields  second  law  using  a  finite  element  technique.  Using  parame¬ 
ter  values  estimated  with  MATLAB  in  the  COMSOL  model,  the  cell 
potential  vs.  DoD  results  were  compared  to  determine  the  error 
between  the  two  programs.  The  comparison  indicated  the  two 
models  differed  by  0.05%  at  most.  Having  established  confidence 
in  the  model  reporting  nearly  identical  results  in  two  solution  pro¬ 
grams,  the  model  is  compared  versus  the  Gomadam  et  al.  model  [1  ] 
to  determine  if  the  new  model  more  accurately  fits  experimental 
data.  Fig.  1 4  plots  cell  potential  versus  time  for  the  1 000  puA  experi¬ 
mental  discharge  data  used  by  Gomadam  et  al.,  the  fit  they  achieved 
to  the  data,  and  the  fit  achieved  with  the  new  model  using  param¬ 
eter  values  in  Tables  4,  5  and  6  .  The  data  for  the  Gomadam  et  al. 
fit  was  taken  directly  from  their  plot  as  the  fit  given  in  their  paper 
was  unable  to  be  replicated  by  the  authors  of  this  paper.  The  new 
model  appears  to  fit  accurately  the  Gomadam  et  al.  discharge  data 
better  than  or  at  least  as  well  as  the  model  reported  in  their  work. 


9716 


D.A.  Strange  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  9708-9718 


Table  4 

Modeling  constants  for  Gomadam  et  al.  [1  ]  cell. 


Parameter 

Units 

Value 

Cli,n,max 

mol  Li+  cur3 

4.825E-02 

Ce 

mol  Li+  cm-3 

IE— 03 

lapp 

A 

-1000E-06 

L 

cm 

5E-05 

ma 

g 

5.109 

M 

- 

15 

Msv 

gmol-1 

595.5 

N 

- 

50 

ne 

mole-  mol-1  SVO 

6 

Sea 

cm2 

2.134E+04 

T 

K 

310.15 

Tref 

K 

310.15 

Pi 

- 

0.5 

P 

gem-3 

4.789 

SLi 

mol  Li+  mol-1  SVO 

6 

Table  5 

Redlich-Kister  coefficients  for  Gomadam  et  al.  [1  ]  SVO  determined  through  param¬ 
eter  estimation.  UfJ  is  in  V  and  the  R-K  coefficient  terms  are  in  units  J  mol-1 . 

Parameter 

Value 

U°„ 

2.902 

Ao 

-2.908E+04 

Ai 

7.992E+03 

A2 

4.055E+04 

A3 

-1.039E+05 

a4 

-1.930E+05 

A5 

7.446E+05 

A6 

3.898E+05 

Ai 

— 2.998E+06 

A% 

1.267E+05 

a9 

6.887E+06 

A10 

— 1.912E+06 

An 

-9.024E+06 

A12 

2.726E+06 

A13 

6.244E+06 

A14 

-1.203E+06 

A15 

-1.733E+06 

Table  6 

Estimated  discharge  fitting  parameters  for  Gomadam  et  al.  [1]  1000  pA  discharge 
data. 


Parameter 

Units 

Value 

cLi,0 

mol  Li+  cm-3 

3.749E-06 

K,  1 

cm2-5  mol-0  5  s-1 

2.081E-09 

kn,2 

cm2-5  mol-0-5  s-1 

9.164E-11 

kn,  3 

cm2-5  mol-0-5  s-1 

9.418E-11 

f\i,l 

cm2  s-1 

7.600E— 15 

f\i,2 

cm2  s-1 

1.405E-15 

Dl  i,3 

cm2  s-1 

5.158E-15 

Rohm,  1 

Q 

13.07 

Rohm,  2 

Q 

10.07 

Rohm, 3 

Q 

7.063 

/Cl 

- 

30 

I<2 

- 

150 

I<3 

- 

100 

I<4 

- 

20 

I<5 

- 

100 

I<6 

- 

30 

#Li,l 

- 

0.100 

#Li,2 

- 

0.340 

@Li,3 

- 

0.300 

#Li,4 

- 

0.540 

#Li,5 

- 

0.150 

#Li,6 

- 

0.300 

A  final  point  of  interest  is  the  effect  of  mechanical  processing 
of  the  SVO  material  on  the  discharge  profile.  This  was  done  for  the 
studied  coin  cell  by  altering  the  diffusion  length  of  the  particle,  L, 
with  results  plotted  in  Fig.  15.  Smaller  particle  sizes  improve  per¬ 
formance  of  the  cell,  but  also  of  note  is  that  particle  length  affects 


Fig.  15.  Effect  of  particle  diffusion  length  on  discharge  profile  at  C/100  (253  pA) 
discharge. 

discharge  behavior  in  the  first  40%  of  discharge.  Were  the  length 
of  the  particle  to  change  during  discharge,  accounting  for  how  the 
length  changes  could  play  a  significant  role  in  modeling.  A  distri¬ 
bution  of  particle  sizes  might  also  exist  as  opposed  to  the  present 
model  which  assumes  all  particles  have  the  same  size. 

4.  Conclusions 

Using  the  R-K  equation  an  excellent  fit  to  experimental  OCP  data 
from  SJM  Li-SVO  cells  was  obtained,  eliminating  the  need  to  use 
empirical  equations  to  model  the  OCP.  Use  of  the  R-K  equation  is 
supported  by  activity  results  consistent  with  literature  sources  and 
industry  standards.  The  single  particle  model  using  Fields  second 
law  can  be  used  to  determine  cell  parameters  (/<n> i ,  DLi ^ ,  DLi  2,  etc.) 
and  also  predict  cell  behavior  during  discharge. 

5.  Future  work 

As  previously  mentioned,  the  method  of  determining  the  dif¬ 
fusion  coefficient  requires  further  refinement.  Examination  of 
discharge  behavior  on  a  different  path  of  intercalation  (i.e.  from 
the  sides)  was  also  considered,  but  was  not  done  in  this  paper. 
While  there  is  evidence  to  support  the  limiting  case  where  lithium 
ions  intercalate  only  through  the  ends  of  the  particle,  no  absolute 
determination  has  been  made.  Examining  the  particle  assuming  a 
limiting  case  where  lithium  intercalates  only  through  the  sides  of 
the  particle  or  relaxing  the  limits  and  letting  lithium  intercalate 
from  all  directions  into  the  particle  may  allow  for  constructing  a 
more  accurate  model. 

Acknowledgements 

The  authors  would  like  to  thank  St.  Jude  Medical  for  supplied 
financial  support  and  experimental  data.  SEM  work  done  at  the  St. 
Jude  Medical  plant  in  Liberty,  South  Carolina. 

Appendix  A. 

The  generalized  electrochemical  reactions  for  intercalation  and 
deintercalation  at  a  solid/solution  phase  interface  as  described  by 
Guo  et  al.  [7]  are  written  as 

T .  ^  DeintercalationT  .+  ^ 

Li  -0S  — >  Li+  +  e-  +  0S 
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Li  -  Os  < —  Li+  +  e  +  Os 

Intercalation 


where  U-0S  is  an  intercalation  site  filled  by  a  lithium  ion  and  0S  rep¬ 
resents  an  unfilled  intercalation  site  [7].  The  concentration  of  filled 
sites  on  the  surface  of  the  particle  is  the  same  as  the  concentration 
of  lithium  ions  on  the  surface,  giving  Eq.  (A.l). 

^U-6S  =  ^Li,j|z=£  —  ^Li,j,sur/^Li,j,max  (A.1) 

Here,  subscript  j  refers  to  the  positive  (p)  or  negative  (n)  electrode. 
The  concentration  of  unfilled  sites  is  then  described  by  Eq.  (A.2). 

—  ^Li,j,max  —  Cu,j|z_£  =  (1  —  xLi ,  j,  surf  )^Li  ,j,  max  (A.2) 

Taking  into  account  the  first  assumption  made  for  the  SP 
approach,  the  electrolyte  lithium  ion  concentration,  ce,  is  also  kept 
constant  at  IE-03  mol  Li+ cm-3  [7]. 

The  current  density  for  the  deintercalation  and  intercalation 
reactions  can  then  be  described  by  Eqs.  (A.3)  and  (A.4),  assuming  a 
negligible  electrolyte  potential. 


h  —  c  ^Lij,max*LiJ,sur/^Xp 

^Li 


Otn  jP  . 


-(01  j  -  02, j) 


(A.3) 


ij  —  kj  c  ^Li  j  maxCe(l  *LiJ,sur/)^xp 
^Li 


®c,jF 


(01 J  -  02 j) 


(A.4) 


where  the  kj  terms  are  the  deintercalation  (-►)  and  intercalation 
(«-)  rate  constants  for  electrode  j,  respectively;  aaj  is  the  anodic 
transfer  coefficient  for  the  deintercalation  reaction  at  electrode  j; 
and  acj  is  the  cathodic  transfer  coefficient  for  the  intercalation 
reaction  at  electrode  j.  An  exchange  current  density  can  also  be 
described  as  the  reaction  current  density  at  open-circuit  conditions, 
or  when  <p,j  -  (p2j  =  Uj. 


I'oj.l  —  max-^Li  j,  sur/ ^Xp 

^Li 


a  a  jF 

a’J  l  J. 

RgT  1 


i(),j,2  —  kj  _  FCijj  maxCe(l  ^LiJ,sur/)^xP 
^Li 

l’oj,i  =  ioj,  2  =  ioj 


ac  iF 
CJ  r  j. 

RgTUj 


(A.5) 

(A.6) 

(A.7) 


If  one  assumes  aaj  +  acj  =  1,  the  exponential  terms  in  Eqs.  (A.5) 
and  (A.6)  can  be  removed  by  raising  Eq.  (A.5)  to  the  acj,  raising  Eq. 
(A.6)  to  the  aaj,  and  multiplying  the  two  together.  If  the  condition 
in  Eq.  (A.7)  is  also  taken  into  account,  Eq.  (A.8)  is  formed. 


ioj  =  =  k“cJk“aj^FcLij-  maxc?aJx%J, 

^Li 


x  (1  Xli, j, surf)  a,J£Xp 


laJuc,j  JT 

RST  J 


exp 


(A.8) 


The  electrode  overpotential,  or  drop  from  open-circuit  condi¬ 
tions  due  to  losses,  is  described  as: 


Oj  =  0i, j  —  02, j  —  Uj  (A.l  2) 

Eq.  (A.l 2)  can  then  be  substituted  back  into  Eqs.  (A.10) 
and  (A.ll).  The  net  current  produced  by  the  intercala¬ 
tion/deintercalation  reaction  is  therefore: 


ij  =  ij  +  ij  =  ioj 


RgT 


exP  l  )  -exP(  --RgfVj 


—  Sl_  ^CLi ,J, max cea’^ surf  (  ^  xLi, j,surf  )  a,j 
oca  jF  \  (  aciF 

6XPI  R^^j-eXP(-R^^ 


(A.13) 


Eq.  (A.13)  is  also  known  as  the  Butler-Volmer  equation  [7]. 

For  the  Li-SVO  reaction  under  consideration  and  neglecting  the 
electrolyte  phase  potential,  the  electrode  overpotential  is  defined 
by  Eq.  (A.14). 


Oj  ~  Ej  Uj 


(A.14) 


Here,  Ej  is  the  electrode  potential  (V).  Because  a  Li/Li+  reference 
electrode  is  used  for  Uj,  the  anode  is  lithium  metal,  and  the  current 
used  is  so  small,  the  anode  potential  is  assumed  to  be  negligible. 
Thus,  the  cathode  potential  plus  the  ohmic  resistance  term  is  the 
cell  potential.  The  anodic  and  cathodic  transfer  coefficients  are  also 
redefined  as  the  symmetry  coefficient,  fy,  multiplied  by  ne/sLi.  With 
the  symmetry  coefficient  held  constant  at  0.5  and  ne/sLi  held  con¬ 
stant  at  one,  Eq.  (A.13)  becomes  Eq.  (A.l 5). 


in  —  fcnFCn  n  maxCe'  XLi, n,sur/¥  xLi, n,surf) 


0.5 


x 


exp 

-exp 


0.5F 


RgT 


( En  ~  Un) 
0.5  F 


RgT 


( En-Un ) 


(A.15) 


The  current  density  is  defined  through  Eq.  (A.16)  as  the  applied 
current  divided  by  the  electroactive  surface  area. 


In  = 


It 


app 

Sea 


(A.16) 


Eq.  (A.15)  can  be  rewritten  as  Eq.  (A.l 7). 


RgT  . 

En  —  lappRohm  +  Un  +  q  ^  0  Slim 


2  /„ 


SE,Fk„ce»-5[cu.0(C(l ,  r)  + 1  )]oa[cL1,„,max  -  cUl0(C(l ,  r)  + 


I)]05}  (A‘ 


(A.17) 


The  exchange  current  density  can  then  be  rewritten  as  Eq.  (A.9), 
assuming  the  rate  constant  kj  =  k^Jk^J  [7]. 

k),j  =  kj—FCi[jmaxEea^X^fjsurjr(j  —  *Li ,j,surf)  a,J  (A.9) 

^Li 

The  current  density  equations  can  be  rewritten  in  terms  of  the 
exchange  current  density. 


ij  =  ioj  exp 


®a,jT 

¥ 


(01 J  ~  02 J  ~  Uj) 


ij  =  -lojexp 


acjF  , 

RgT 


(01  j  -  02  j  -  Uj) 


(A.10) 

(A.11) 


Appendix  B. 

An  examination  of  electrochemical  potentials  allows  for  an 
alternative  method  for  estimating  the  OCP  of  lithium  intercala¬ 
tion  electrodes,  though  in-depth  analyses  of  these  potentials  are 
beyond  the  scope  of  this  work.  The  method  replicated  here  was 
developed  by  Karthikeyan  et  al.  [13],  restricted  in  scope  in  this 
paper  to  estimating  the  OCP.  First,  the  intercalation  electrode  equi¬ 
librium  potential  relative  to  a  lithium  reference  electrode  can  be 
represented  by  Eq.  (B.l),  where  a  refers  to  the  intercalated  species 
and  p  refers  to  vacant  sites  [13]. 

—FUn  =  /Xy  +  Up  -  [Tol 
*u 


(B.l) 
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Here,  the  electrochemical  potential  of  species  i  can  be  expressed 
in  terms  of  an  activity  coefficient,  mole  fraction,  and  reference  elec¬ 
trochemical  potential  via  Eq.  (B.2). 

Mi  =  if  +  RgT  ln(y,x,)  (B.2) 

where  subscript  i  refers  to  a  or  (3.  Eq.  (B.2)  can  be  substituted  into 
Eq.  (B.l)  resulting  in  Eq.  (B.3),  which  can  be  used  to  represent  the 
OCP  profile. 

-  <4  +<4 -< tvm  ©  +W"  © 

.|FU«  +  R„n„g)+Ejn„(a)  (B.3) 

If  unit  activity  is  assumed,  Eq.  (B.3)  reduces  to  the  classic  Nernst 
equation,  Eq.  (B.4)  [13]. 

~~FUn  =  ^FU°  +  RgTln  (B.4) 

sLi  sLi  V  Xa  J 

If  unit  activity  is  not  assumed,  the  activity  coefficients  can  be 
related  back  to  the  partial  molar  excess  Gibbs  free  energy  through 
Eq.  (B.5). 

RgTln(yi)  =ff  =  ^-(ntCE)TA,  +  {  (B.5) 

Here,  subscript  l  refers  to  a  or  (3.  The  R-K  equation  can  then  be 
used  to  approximate  the  excess  Gibbs  free  energy  [13]. 

M 

Gf  =  Xq^nXpjt^  ^k(Xq,n  ~  X^  n) 
k= 0 


M 


—  Xoi,n{^  ~  Xpi^n)^  ^-Afc(2Xo;,n  —  1) 


(B.6) 


k= 0 


Eq.  (B.6)  can  be  substituted  back  into  Eqs.  (B.3)  and  (B.5)  to 
obtain  Eq.  (B.7).  Eq.  (B.7)  is  the  form  used  in  Eq.  (21 )  for  the  OCP. 

fwn  =  fFU°n  +  RgTref  In  ( 1~XLi'n’su,f 

^Li  ^Li  y  Xi[n,surf 

J  rnv  i  -ik+1  2xLi  „  surf k(l  -XLi  n  SUrf) 

+  <  ^JKk  •  (2xLi,n,sur/  -  1) - — -  .i_k 

\<= o  L  l^Li ,n,surf  U 

(B.7) 


One  can  also  solve  for  the  activity  coefficient  of  the  intercalated 
material  via  Eq.  (B.8)  [13]. 

M 

RgT  ln(yu,n,  surf)  —  ^  }^/c(l  —  ^Li,  n,surf)  (2^Li,  n,surf  ~ 
k= 0 


x 


2x. 


Li,  n,  surf  * 


(2xLi ,n,sur/  -  1) 


+  1 


(B.8) 


Appendix  C. 


The  treatment  of  the  governing  equation,  Eq.  (4),  and  boundary 
conditions,  Eqs.  (6)  and  (7),  while  utilizing  MATLAB  deserves  men¬ 
tion,  as  MATLAB  does  not  automatically  solve  for  Fields  second  law 
whereas  COMSOL  does.  Instead,  Fields  second  law  is  broken  apart 
from  a  partial  differential  equation  into  a  series  of  ordinary  differ¬ 
ential  equations  using  the  method  of  lines  approach.  This  divides 
the  particle  along  its  length  into  individual  nodes,  calculating  the 
concentration  at  each  point.  For  the  modeling  performed  in  this 
work,  N  =  50  nodes  are  used.  Defining  the  center  of  the  particle  as 
the  first  node,  m  =  1,  and  the  surface  of  the  particle  as  the  last  node, 
m  =  N=  50,  Eq.  (C.l)  solves  for  the  central  nodes. 


dCti  ,n,m 

dt 


=  DLi 


(^Li,n,m-1  n,m  +  ^Li,n,m+1 ) 

(L/(JV-1))2 


m  =  2,  N  -  1  (C.l) 


where  L/(N  -  1 )  breaks  the  particle  into  equal-length  segments  and 
L  is  the  half-length  of  the  particle  (cm).  The  boundary  conditions 
at  the  center  of  the  particle,  Eq.  (C.2),  and  at  the  surface,  Eq.  (C.3), 
are  also  included,  using  forwards  and  backwards  finite  difference 
techniques,  respectively. 


n  \  ^CLi, n,N 

^Li(cLi  ,n,iv) — ^ — 


^Li(cLi  ,n,iv) 


3Cu,n,N  ~  4cli, n,N- 1  +  cLi,n,N-2 

2(L/(N  —  1)) 


Iqppsu  =  q 

neFSEA 


(C.3) 
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